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Abstract 

Consider “Frozen Random Walk” on Z: n particles start at the origin. At any 
discrete time, the leftmost and rightmost [j J particles are “frozen” and do not move. 
The rest of the particles in the “bulk” independently jump to the left and right uniformly. 
The goal of this note is to understand the limit of this process under scaling of mass 
and time. To this end we study the following deterministic mass splitting process: 
start with mass 1 at the origin. At each step the extreme quarter mass on each side 
is “frozen”. The remaining “free” mass in the center evolves according to the discrete 
heat equation. We establish diffusive behavior of this mass evolution and identify the 
scaling limit under the assumption of its existence. It is natural to expect the limit 
to be a truncated Gaussian. A naive guess for the truncation point might be the 1/4 
quantile points on either side of the origin. We show that this is not the case and it is 
in fact determined by the evolution of the second moment of the mass distribution. 


1 Introduction 

The goal of this note is to understand the long term behavior of the mass evolution process 
which is a divisible version of the particle system “Frozen Random Walk”. We define 
Frozen-Boundary Diffusion with parameter a (or FBD-q) as follows. Informally it is 
a sequence fit of symmetric probability distributions on Z. The sequence has the following 
recursive definition: given fit., the leftmost and rightmost ^ masses are constrained to not 
move, and the remaining 1 — a mass diffuses according to one step of the discrete heat 
equation to yield fit+i- In other words, we split the mass at site x equally to its two 
neighbors. Formal descriptions appear later. We briefly remark that this process is similar 
to Stefan type problems, which have been studied for example in |2j. 

Now we also introduce the random counterpart of FBD-a. We define the frozen random 
walk process (Frozen Random Walk-(n, 1/2)) as follows: n particles start at the origin. At 
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any discrete time the leftmost and rightmost particles are “frozen” and do not move. 
The remaining n — particles independently jump to the left and right uniformly. 

Letting n —> oo and fixing f, the mass distribution for the above random process converges 
to the t th element, /x*, in FBD-a. However, if t and n simultaneously go to oo, one has to 
control the fluctuations to be able to prove any limiting statement. Figure [l] depicts the 
mass distribution nt and the frozen random walk process for a = ^. 



Figure 1: Frozen-Boundary Diffusion-^ and Frozen Random Walk-(10000, aver¬ 
aged over 15 trials at t = 25000. 


At every step t of FBD-ct, we also keep track of the location of the boundary of the 
process, fit, which we define as 

Pt := sup jx E Z : m ([x,oo)) > ^ j . 

We will show that 

Lemma 1. For every a E (0,1) there exist constants a, b > 0 such that 

ay/t < j3t < b\ft , Vi. 
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The lemma above suggests that a proper scaling of fit is \ft. Motivated by this behavior 
of the boundary fit, one can ask the following natural questions: 

Ql. Does ^j= t converge? 

Considering y t as a measure on M, for t = 0,1,... define the Borel measure fifiot) = jit 
on M equipped with the Borel a— algebra such that for any Borel set A. 

P-t(A) = fi t {{yVt : y e A}). (1) 

We can now ask 


Q2. Does the sequence of probability measures yt have a weak limit? 
Q3. If yt has a weak limit, what is this limiting distribution? 

We conjecture affirmative answers to Ql and Q2: 

Conjecture 1. For every a G (0,1), there exists £ a > 0 such that 

lim = i a . 

t-too y/t 


Conjecture 2. Fix a G (0,1). Then there exists a probability measure Moo(c>0 on M such 
that as t —>• oo, 

~ weak , s. 

yt = =>■ Moo (a), 

where denotes weak convergence in the space of finite measures on M. 

That Conj. [2] implies Conj. [T] is the content of Lemma |4j We now state our main result 
which shows that Conj. [l] implies Conj. [2] and identifies the limiting distribution, thus 
answering Q3. To this end we need the following definition. 


Definition 1. Let <&(•) be the standard Gaussian measure on M. Also for any q > 0 denote 
by d)q(-), the probability measure on M which is supported on [—q.q] and whose density is 
the standard Gaussian density restricted on the interval [—q, q] and properly normalized to 
have integral 1. 


Theorem 2. 


Assuming that lim -^= 

t-t OO y/t 


is a constant, the following is true: 


~ weak f n 
fit -^ Moo(^)j 


where, 

Oi Of. 

Moo (a) = -^S(-qa) + (1 - Oi)^ qa + -5(q a ), 
and q a is the unique positive number such that: 

a (1 — 

2 Qa VTk $([-&*,&*]) 

Remark 1. It is easy to show (see Lemma [4]) that the above result implies that 

Pt 

lim —n= ( lo l - 
t-t oo y/t 

Thus observe that by the above result, just assuming that the boundary location properly 
scaled converges to a constant determines the value of the constant. This is a consequence 
of uniqueness of the root of a certain functional equation discussed in detail in Section [2j 
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1.1 Formal definitions 


Let FBD-a := {/zq, /zi ,...}: where for each i = 0,1, ...,//* is a probability distribution on 
Z. For brevity we suppress the dependence on a in the notation since there is no scope of 
confusion as a will remain fixed throughout any argument. 

Let /ro = 5(0) be the delta function at 0. By construction yt will be symmetric for all 
t. As described in Section [I] each y t contains a “constrained/frozen” part and a “free” 
part. Let the free mass and the frozen mass be denoted by the mass distributions vt and ft 
respectively. 

Recall the boundary of the process, 


Then for all y > 0, 


Pt = sup jz G Z : fjL t ([z,cx>)) > ^ j 


ft(y) ■■= 


My) y> Pt 
y = Pt 

z>Pt 

0 otherwise. 


( 2 ) 

(3) 


For y < 0 let ft{y) ■= ft(~y )• Thus ft is the extreme a/2 mass on both sides of the origin. 
Define the free mass to be vt '■= yt — ft- With the above notation the heat diffusion is 
described by 

v t [x- l) + v t (x + l) 

yt+ i(s) =-g- + ( 4 ) 

Recall Lemma [TJ which implies the diffusive nature of the boundary: 


Lemma 1. For every a G (0,1) there exist constants a,b > 0 such that 

ay/t < Pt < by/t ,Vi. 

This result implies that in order to obtain any limiting statement about the measures fit, 
one has to scale space down by y/t. 

The proof of the lemma appears later. Let us first prove that the frozen mass ft cannot 
be supported over many points. 

Lemma 3. For all t, the frozen mass at time t, ft, is supported on at most two points on 
each side of the origin, i.e., for all y G Z such that \y\ > Pt + 2, we have ft(y) = 0. 

Proof. The lemma follows by induction. Assume for all k <t, for all y such that \y\ > Pk+ 2, 
we have yk{y) = fk{y) = 0. The base case t = 0 is easy to check. Now observe that by Q 
and the above induction hypothesis, 

ht+i{y) = 0, (5) 

for all \y\ > Pt + 2. Also notice that by Q it easily follows that Pt is a non-decreasing 
function of t. Thus clearly for all y, with |y| > p t +1 + 2 > p t + 2, 

y t +i(y) = o. 


Hence we are done by induction. 

O' 
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We now return to the proof of the diffusive nature of the boundary of the process fit- 


Proof of Lemma [l} We consider the second moment of the mass distribution /// . which 
we denote as M 2 (t) := ^^/it(x)x 2 . This is at most (fit + l) 2 since /it is supported on 

X€Jj 

[—/3t — 1, fit + 1] by Lemma [ 3 J It is also at least a fit 2 since there exists mass a which is at a 
distance at least fit from the origin. Now we observe how the second moment of the mass 
distribution evolves over time. Suppose a free mass m at x splits and moves to x — 1 and 
x + 1. Then the increase in the second moment is 

^((x + l) 2 + (x — l) 2 ) — mx 2 = m. 

Since at every time step exactly l —a mass is moving, the net change in the second moment 
at every step is 1 — a. So at time t the second moment is exactly 


t( 1 — a). 


( 6 ) 


Hence a fit 2 < t(l — a) < (fit. + l) 2 , and we are done. E[ 

We next prove Conjecture [T] (a stronger version of Lemma [I]) assuming Conjecture [2j 


Lemma 4. If Conjecture [I| holds, then so does Conjecture [TJ i.e., 
exists l a > 0, such that 


lim = l D 

i-s-oo y/1 


for every a G (0,1), 


there 


Proof. Fix a G (0,1). From Lemma [I] we know that {fit./y/t} is bounded. Hence, if fit/Vt 
does not converge, there exists two subsequences {si, S 2 , ■ ■ •} and {t\,t 2 , ■ ■ ■} such that 


lim fisi/y/si = qi and lim fitjy/tj -t 

i— > 00 j— ^00 

for some q 2 . q\ > 0 such that q 2 — q\ '.= 5 > 0. Recall /ioo (a) from Conjecture [2| Now by 
hypothesis, 

lim n 8i = ii oo(«), I™ Hu => HooKpL). 

i—t 00 

This yields a contradiction since the first relation implies /ioo (ct) assigns mass 0 to the 
interval (92 — f, q 2 + f) while the second one implies (by Lemma [ 3 ]) that it assigns mass at 
least j to that interval. O 


2 Proof of Theorem [2] 

The proof follows by observing the moment evolutions of the mass distributions /it and using 
the moment method. The proof is split into several lemmas. Also for notational simplicity 
we will drop the dependence on a and denote /ioo (a) and q a by /i and q respectively since 
a will stay fixed in any argument. Thus 

M := /Xoo(a) = f S(-q) + (1 - «)$, + ^6(q). (7) 

Also denote the k th moments of /it as M^(t). We now make some simple observations 
which are consequences of the previously stated lemmas. Recall the free and frozen mass 
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distributions Uf and ft . We denote the k th moments of the measures vt (the free mass at 
time t), ft (the frozen mass at time A), by M'f (A) and M k (t) respectively. Also define ft and 
i*t similarly to fit in ([!]). Assuming Conjecture [TJ it follows from Lemma [3] that, 


x weak 

Jt 


f 


( 8 ) 


where / := fS(—£ a ) + ^S(£ a ), and £ a appears in the statement of Conj[l| This implies that 

M f k (t) _ f 0, k odd 

t k / 2 \ a£^,(l + o(l)), A: even 


where o(l) goes to 0 as t goes to infinity. 


The proof of Theorem [2] is in two steps: first we show that £ a = q and then show that 
vt converges weakly to the part of /r which is absolutely continuous with respect to the 
Lebesgue measure. Clearly the above two results combined imply Theorem [2j 


As mentioned this is done by observing the moment sequence Now notice owing to 

symmetry of the measures for any t , M^k+iit) = 0 for all non-negative integers k. 

Thus it suffices to consider M 2 k(t ) for some non-negative integer k. We begin by observing 
that at any time t the change in the moment M- 2 k(t + l) — M 2 k(t), is caused by the movement 
of the free mass i't- The change caused by a mass m moving at a site x (already argued in 
the proof of Lemma [l] for k = 1) is 


m((x + l) 2fc + (x 
2 



( 10 ) 


Now summing over x we get that, 


k , 

M 2k (t + 1) - M 2k (t) = ^2 

i=i ^ 


2k 

2k — 2 i 


^2k-2i^)- 


(ID 


Notice that the moments of the free mass distribution ut appear on the RHS since m in (10) 


was the free mass at a site x. Now using we sum M 2 k(j + 1) - M 2k (j) over 0 < j < t- 1 
and normalize by t k to get 


M 2 k(t) 
t k 



( 12 ) 


Recall that by LemmajlJ for any k > 1, M!f k _ 2 (jf) is 0(j k x ). Moreover, the above equation 
allows us to make the following observation: 

the existence of lim 

00 


Claim. Assume ([9]) holds. Then for any A: > 1, 


M 2k -2^i) 

Ak- 1 


implies 


existence of lim 

T+oo 


M&U) 

jk 


Proof of claim. Notice that by Lemma 0 Mfk_({j) = 0(j k 2 ) for any £ < 4. Also let 


ML Jj) 

T Zk — 2\ J / 71 /rV 

lim -^3]— = M 2k _ 2 , 

J— ^OO J ^ L 
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which exists by hypothesis. 

Thus using (12) and the standard fact that 


l _1 Ak—X i r 1 

lim -j-r- = / x^ 1 

t ^ oo ~' 0 t f Jo 


dx = — 
k 


we get 


Thus 


t-1 

E 


i=i Li=i 


E 


\ 


2k-2iJ t k 

M2k(t) 


— (2k — l)M 2k _ 2 + o(l) + 0 ( - ) . 


lim 

t— >oo 


and since 

we are done by ([9]). 


^ = (2fc-l)M&_ 2 

M 2fc (t) = M&(t) + M/ fc (f), 


(13) 

(14) 

□ 


Using the above claim, the fact that lim ^ and hence, lim k 1 ^ (by {fl4|) and d9l)) 

t —^oo t k / 2 t->oo t k / 2 - ^ 

exists for all /c, follows from the fact that M |^ = (1 — a) (see ©»■ Let us call the limits 

M k and respectively. 


Thus we have 

M 2fc = M 2k + = (2k - 1 )M 2 V 2j (15) 

where the hrst equality is by (14) and © and the second by ( |13[ ). For k = 1 we get 

al 2 0 + M 2 = 1 - «• 

Notice that this implies that for all k, ML, can be expressed in terms of a polynomial in 
l a of degree 2k, which we denote as P k (l a )- Then, by (15) the polynomials P k satisfy the 
following recurrence relation: 


P k (l a ) = (2k~l)P k _ 1 (lo)-al 2k 
P 0 = 1 — a. 


(16) 


By definition, we have 


Pkita) = M 2k = lim = lim 

t—> OO t 


t—¥ OO 


E 

-pt<x<pt 


(17) 


Thus assuming Conj. 
family of inequalities 


and the fact that z^(x) = 1 — a for all t, we get the following 


0<P k (l a )<(l~a)l 2 a k Vk>0. (18) 

We next show that the above inequalities are true only if l a = q where q appears in Q . 
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Lemma 5. The inequalities in (IS) are satisfied by the unique number £ a such that 

a (1 — a)e -f “/ 2 

2 ““ v^([- 4 , 4 ]) 

where <&(■) is the standard Gaussian measure. 

Thus the above implies that necessarily £ a = q where q appears in Q. This was men¬ 
tioned in Remark [l] 

Proof. To prove this, first we write the inequalities in a different form so that the polynomials 
stabilize. To this goal, let us define 


Pi, = 


Pk 


{2k — 1)!! 


where (2 k — 1)!! = (2 k — 1)(2 k — 3)... 1. Then it follows from (16) that 


h (4) = Pk- i(4) - 


a 


(2 k - 1)!! 


(/2k 

-a • 


Hence 


The inequalities in (18) translate to 


n ,. n y < , e +2 

{2i - 1)!! “ {2k + 1)!! 


(19) 


Let us first identify the power series 

9{x) = 


00 x 2i_1 


i =1 


(2 i - 1)!!' 


Clearly the power series converges absolutely for all values of x. It is also standard to show 
that one can interchange differentiation and the sum in the expression for g{-). Thus we get 
that, 

dg{x) i , /x 

Solving this differential equation using integrating factor e~ x “/ 2 and the fact that g{ 0) = 0 
we get 


g{x) = e 




f 


-y 2 / 2 dy. 


As k —> oo, the upper bound in (19) converges to 0 for any value of £ a . Also the expression 
in the middle converges to 1 — a — cd a g{£ a ). Thus taking the limit in (19) as k —> oo we 
get that £ a > 0 satisfies 

e a g(£ a ) = — ■ (20) 

a 

Clearly this is the same as the equation appearing in the statement of the lemma. Also 
notice that since xg{x) is monotone on the positive real axis, by the uniqueness of the 
solution of (20) we get £ a = q where q appears in Q. Hence we are done. D- 
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The value of £ a that solves (20) when a = 1/2 is approximately 0.878. Figure [2] shows 
the numerical convergence of fit./Vt to q a for various values of a. 


% 

cqT 


1.736 

1.193 

0.878 

0.615 


1000 2000 3000 

t 


a=0.7 

a=0.5 


a=0.3 

a=0.1 


4000 


5000 


Figure 2: Convergence of fit/Vt for various a. The horizontal lines denote the values q a 
and the curves plot as a function of time t. 


Thus, assuming Conjecture [lj by Lemma [HJ ft converges to / (as stated in Q) which 
consists of two atoms of size ^ at q and —q. To conclude the proof of Theorem [2j we now 
show v t converges to the absolutely continuous part of /i (see 0). Recall that by © and 
Lemma [ 5 ] the 2k th moment of i\ converges to Pk(q). We will use the following well known 
result: 

Lemma 6 (30.1, [I]). Let /r be a probability measure on the line having finite moments 
Oik = f- ^ x k p(dx) of all orders. If the power series akr k /k\ has a positive radius of 

k 

convergence, then fa is the only probability measure with the moments a\,a 2 , ■ ■ ■■ 

Thus, to complete the proof of Theorem [2] we need to show the following: 

Claim. The 2 k th moment of the measure (1 — a)& q is Pk{q) where q appears in Q. 


To prove this claim, it suffices to show that the moments of (1 — a)& q satisfy the recursion 
(16). Recall that q = q a . Let C = C a := ^ Using integration by parts we have: 


x"e 


1 „2k n 2 d x 


C 


f 9 

J-Q 




xe 2 dx 


C 


2 q 


,2k lg 2 


c 


+ 


(U-D l\^ e -4 ix . 


c 


<-q 
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By the relation that q satisfies in the statement of Theorem [2j the first term on the RHS 
without the — sign is aq 2k . Also, note that the second term is (2 k — 1) times the (2k — 2) nd 
moment of (1 — a)Q q . Thus, the moments of (1 — a)Q q satisfy the same recursion as in (16). 


Now from Example 30.1 in [I], we know that the absolute value of the k th moment of 
the standard normal distribution is bounded by k\. Then, similarly, the absolute value of 
the k th moment of our truncated Gaussian, & q , is bounded by c k k\ for a constant c. Then 
Lemma [6] implies that is determined by its moments and quoting Theorem 30.2 in [1J 
we are done. 

□ 


3 Concluding Remarks 

We conclude with a brief discussion about a possible approach towards proving Conjectures 
[l] and [2] and some experiments in higher dimensions. 
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x 


Figure 3: Heat map of the free mass distribution after 1000 steps in 2 dimensions for 

FBD-1/2. 

The free part Vt. of the distribution fit could represent the distribution of a random walk 
in a growing interval. If the interval boundaries grow diffusively, the scaling limit of this 
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random process will be a reflected Ornstein-Uhlenbeck process on this interval [—q,q]. We 
remark that the stationary measure for Ornstein-Uhlenbeck process reflected on the interval 
is known to be the same truncated Gaussian which appears in Theorem[2j see [3] (31)]. This 
connection could be useful in proving the conjectures. 

We also note that similar results are expected in higher dimensions; in particular, the mass 
distribution should exhibit rotational symmetry. See Fig [3} Note the truncated Gaussian 
shape for the slices x = 0 (Fig@. 



Figure 4: Slice of the free mass distribution at x = 0 after 1000 steps in 2 dimensions for 
the analogue of FBD-1/2. 
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